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Abstract 

We analyze the effects of tlie environment on the spin tunneling process of param- 
agnetic and superparamagnetic particles and conclude that the "bare" macroscopic 
tunneling rate is hardly affected in such case, but others more effective processes 
come up to change the magnetization state of the particle. We measure the degree of 
coherence loss by the linear entropy which will account for all correlations present in 
the process. We conclude that for both, paramagnetic and superparamagnetic par- 
ticles, the decoherence time scale is extremely short {'^ indicating that 
coherent tunneling should be strongly suppressed in favor of incoherent tunneling, 
i.e., the population of higher levels with subsequent decay. Interestingly enough, 
the ground state tunneling rate is hardly affected by such dissipative mechanisms. 
However other processes immediately give rise to new (incoherent) possibilities of 
crossing the barrier. 
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Recently discovered magnetic molecules[l,2,3,4] are of interest because of the 
gigantic relaxation times of their magnetization. Recent experiments on mag- 
netic relaxation of molecular crystals of Mni2 and Fes find strong evidence for 
tunneling mediated relaxation at low temperatures[5,6,7,8,9,10]. Several theo- 
retical studies have been performed[ll,12,13,14], leading to tunneling rates of 
the order of a few months. The situation with superparamagnetic particles is 
less clear, however. A change in magnetization direction at low temperatures 
is indeed observed, the question remaining as to whether this is due to a uni- 
tary tunneling effect. Contradictory statements regarding the observability of 
the tunneling effect can be found in the literature[15,16]. 
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In the present letter, we set up a collective model for the tunneling of mag- 
netic molecules, derive the corresponding tunneling rates and include dissipa- 
tion via coupling with a Caldeira-Leggett type of environment. As expected, 
the decoherence time scale is found to be much shorter than the characteris- 
tic tunneling time. However, surprisingly enough, we find that the tunneling 
rate is hardly affected by the environment induced dissipation. We therefore 
conclude, within the scope of the present model, that the spin-phonon inter- 
action does not significantly interfere with the coherences which are essential 
for the tunneling process. It however very quickly opens up several alternative 
routes for the system to change its magnetization, rendering the unitary tun- 
neling mechanism but a possibility among many others. Our result strongly 
supports the picture of incoherent tunneling, since coherence is destroyed just 
about immediately. 

We start by considering a Hamiltonian which describes the magnetization 
properties of such molecules [17], 



The anisotropy constants satisfy A ^ B > 0. and S± = iSy are spin 
operators. The energy sphtting of the ground state due to the second term on 
the r.h.s. of Eq. (1) can be calculated by means of path integral techniques 
and the instanton method. Details of this calculation can be found in Ref. [18]. 
We get 



This result can be compared with that obtained by Hartmann-Boutron by 
means of others methods[19]. For magnetic molecules, like Mni2, s = 10 and for 
typical superparamagnetic particles, s ~ 3000/i. For both we have B <^ A, and 
therefore AE'inst ~ 0, suggesting then that th observation of this phenomenon 
in the laboratory is impossible. 

Obviously thermally activated over-barrier turnover of the magnetization di- 
rection at sufficiently high temperatures (normal superparamagnetism) is pos- 
sible. In the following, however, we analyze the possibility of spin-tunneling, 
i.e., under-barrier changes of the magnetization direction. 

The model of a single collective variable is apparently too naive. We therefore 
next take into account the fact that the magnetic particles are not isolated and 
estimate the characteristic time of environmental effects. We consider phonon 
degrees of freedom, represented as a set of harmonic oscillators described by 
the Hamiltonian 



Hs^-ASl + B{Sl + S^_). 
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and coupled to Hs (this model has been proposed in [20]) 



//int = BUSl + SV) E sjl^kM + (4) 

where Sint is an overall interaction constant, M the particle mass, ka the 
wave number associated to the phonon with frequency c<Jq and a]^[aa) creation 
(annihilation) operators of the environmental degrees of freedom. 

Assuming this interaction to be weak enough with respect to the collective 
dynamics, Eq. (1), we use perturbation theory to evaluate the reduced density 
of the spin degree of freedom. 



p(^)=p(o)--y^ dt'[H,^,{t'),~pm 

_J_ fdt' f dt"[^int(OJ^int(O,P(0)]] + 

n JO JO 

= p(0)+^(l)+^(2) + ... 
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where ^int = e'^"* Hi^t e *^o* and Hq = Hs + Hr. The reduced density is then 

Ps = Trij[p]. (6) 

Tr^ represents the trace over the reservoir (phonons) variables. From the above 
expressions it is not difficult to obtain the two quantities which lead us to our 
final conclusions: 

(1) The tunneling probability, i.e., the probability as a function of time to 
find our system in state | — s) , having been initially prepared in | + s) 



p.s{t) = i:v[\-s){-s\ps{t)] 
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(2) The coherence loss as measured by the linear entropy (or idempotency 
defect [21]) 



5s{t) = 1 - l^s[{ps?] 

S™2(^m,m+2)^ 



+ 



(8) 



From this formula we also define the "decoherence time" (rdec) as being 
the inverse of the factor multiplying t, or in others words, the time at 
which 5s{t) — 1. 
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In the above expressions, we have used a spectral function for the phonons 
given by 

and a hnear dispersion relation for the phonons cu^ = cka, where c is the sound 
velocity in the medium. Moreover, Debye's model is used to transform the sum 
into an integral, i.e. 



Z^Qi^a)^ dujg{uj)Q{uj) 



with g{u!) = 9a;^/a;|) and cu^^ = 3V/{18n^c^), V being the particle's volume. 
Nevertheless it is only necessary to know the value of the density p since V 
and M always appear in formulas in the form M/V . With these ingredients, 
Eq. (9) becomes 

corresponding to superohmic dissipation, according to Caldeira and Leggett's 
model[22]. In Eq. (8), s'l^ = (s ^ m)(s ± m + ^ m - ± m + 2) and 
<^m,m±2 = {Ejn - Em±2)/fi. In Eq. (7), S-2 IS cqual s!!'2> for s. 

Let us begin by analyzing the tunneling rate. The environmental effects are 
contained in the brackets of Eq. (7) and are, for typical values of the constants 
(see Table 1), much smaller than 1. We can therefore conceive that the terms 
in brackets correspond to the first terms in the expansion of the exponential 
e-7(/3)t^ with 

This result is thus very analogous to the tunneling rate obtained in the context 
of the spin-boson model, i.e., the "free" tunneling frequency is damped. Our 
model hypothesis, which we believe valid for magnetic particles, indicate that 
the tunneling rate remains essentially unaffected by the presence of the other 
degrees of freedom. This means, in particular, that the unitary quantum me- 
chanical tunneling process will be hardly affected by this interaction. However, 
a careful look at the dccoherence time (see Fig. 1) shows that very quickly 
(Tdec ~ 10~^^s (~ 10~^s) for superparamagnetic (Mni2) particles), many new 
channels will be open, the density matrix will be "contaminated" by other 
available states and lose purity. Physically this means that, although the par- 
ticular coherences necessary for the macroscopic quantum coherent tunneling 
of the ground state remains essentially unaffected by the dissipation, states, 
candidate to tunnel the barrier, will not do it in a coherent way since others 
mechanisms for magnetization change will come into play very quickly, thus 
giving support to the idea that the experimentally observed magnetization 
change must be phonon assisted. 
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Table 1 

Typical parameter values of the Mni2 molecule and a superparamagnetic particle 
Parameters Mni2 molecule^^^ Superparamagnetic particle^^-* 

A 7.5 X lO-^'^J 2.57 x IQ-^^ J 

B 1.7 X ~ IQ-^A 

Bi^t ^A w 13 cm.-^ = 4, 1 X lO^^s j 

s 10 3222 

density 1.83 x 10^ kg/m^ 5.0 x 10'^ kg/m^ 

sound's velocity 2.0 x 10'^ m/s 3.0 x 10'^ m/s 

Sources: Leuenberger and Loss [23]; ^"^^ Pfannes et al.[16] 
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Fig. 1. Inverse of the decoherence time as function of the initial state for the Mni2 
(a) and a typical superparamagnetic particle (b) at different temperatures. 
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